Sterols sense swelling in lipid bilayers 
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In the mimetic membrane system of phosphatidylcholine bilayers, thickening ( "pre-critical behavior", 
"anomalous swelling") of the bilayers is observed, in the vicinity of the main transition, which is non- 
linear with temperature. The sterols cholesterol and androsten are used as sensors in a time-resolved 
simultaneous small- and wide angle x-ray diffraction study to investigate the cause of the thickening. We 
observe precritical behavior in the pure lipid system, as well as with sterol concentrations less than 15%. 
To describe the precritical behavior we introduce a theory of precritical phenomena. The good temperature 
resolution of the data shows that a theory of the influence of fluctuations needs modification. The main 
cause of the critical behavior appears to be a changing hydration of the bilayer. 



I. INTRODUCTION 
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2 Our ultimate aim is a fuller description of biological cell 
Qmembranes, via the physics of membrane models. Cell 
,^rnembranes are mainly in equilibrium, yet obviously must 
^-e xperience fluctuating non-equilibrium conditions during 
Qcellular events such as membrane fusion and transmem- 
Obrane signaling. An excellent experimental model system 
for biological cell membranes is a dispersion of phospho- 
lipid in water, known as liposomes, which exhibit a variety 
;>of interesting phases as temperature is changed |Q. Ma- 
Cjor changes in physical properties are seen at the gel to 
ONluid phase transition, the main transition at temperature 
^iTm, which is a first-order transition P|§, occurring over 
*• narrow temperature range. However, when approach- 
OQing the main transition from above T m , structure-sensitive 
Q^methods reveal a nonlinear increase (with decreasing tem- 
~p<perature) of the lattice parameter d in phosphatidylcholine 
c3(PC)-membranes which is attributed to critical be- 

Qiavior (illustrated in Fig. ||). This increase is here termed 
I "anomalous" behavior. This term has been used also in 
""Ojther contexts || . The corresponding critical temperature 
would lie at a temperature just below that of T m , and 
^o is "hidden" by the main transition. Critical behavior 
• -is also seen in liposomes by nuclear magnetic resonance 
t ^(NMR) J5J and ultrasonic absorption S . This behavior in 
^^phosphatidylcholine membranes, a major component of cell 
{^membranes, is believed to be of physiological importance, 
^because small temperature changes change the thickness of 
the membrane drastically and thus can easily affect protein 
functions Jj|[lO| . Additional support for biological relevance 
comes from the fact that the nonlinear behavior occurs in 
the neighborhood of a phase transition, a non-equilibrium 
state of matter. 

The bilayer thickness d, measured by x-ray diffraction, is 
composed of three regions which are the headgroup thick- 
ness (iff, which is part of the bilayer thickness ds, and the 
inter-bilayer water layer dw, where ds + dw — d. The 
aim is to disentangle the contribution of each region to the 



observed critical thickening in d, by examining four physi- 
cally different models, each of which affect mainly one re- 
gion; the models do not necessarily agree In brief, 
Model I (following the description of the models by Na- 
gle [pi) attributes the thickening to critical fluctuations, 
where a softening of the bilayer leads to a reduced bend- 
ing modulus (via enhanced Helfrich undulations) [ |l2[ , and 
affects mainly dw |l^7@ • This model has been tested and 
supported by neutron scattering analysis ||l1J. An anal- 
ysis of x-ray data Q did not confirm the reduced bending 
modulus. Model II attributes the thickening to an increas- 
ing conformational order in the hydrocarbon chains upon 
approaching T c , and hence to an increase in the lipid bi- 
layer thickness ds 0- Both analyses use mean field- type 
liquid crystal theories whose application to critical phe- 
nomena and phase transitions - according to the Landau 
theory the L^-P^-transition has no chance of being second 
order jl5j - is debatable. In Model III, the thickening is at- 
tributed to interbilayer hydration and van der Waals forces, 
which increase the water layer thickness dw [^)- Reasons 
for a change in these forces when approaching the transition 
and their effect on the bilayer are not well understood, nor 
is the functional form modeling these forces tested in the 
neighborhood of the main transition. Model IV attributes 
the thickening to an increased headgroup layer thickness 
du only, which is, however, not supported by NMR work 
& 

Sterols are an important constituent of eukaryotic mem- 
branes. In mammals, cholesterol is essential for the well- 
being of cells H ; a concentration-dependent effect on soft- 
ening and stiffening the membrane has been found for 
cholesterol Jll| . Sterols can also play a role in some protein 
functions coupling lipid dynamics to protein dynamics Q , 
and in the action of an acetylcholine receptor p [.Choles- 
terol organizes submicron domains in living cells] IS . Here, 
in the absence of proteins, sterols serve us as a natural sen- 
sor of the lipid dynamics. In this context it is worthwhile 
to mention that a high content of cholesterol (> 20 mol%) 
changes the lipid phase diagram drastically, whereas small 
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amounts (< 10 mol%) do not affect the lipid phase dia- 
gram significantly ||| . Note that a phase diagram contains 
information about the equilibrium phases, and that it does 
not give direct information about non-equilibrium effects 
which may naturally occur at the phase boundaries. It has 
been shown |ll[] that small amounts of cholesterol strongly 
affect the structure of the lipid membrane near the main 
transition. 

Our first aim in studying the phenomenon of increased 
swelling is to find out whether it is analytical or indeed 
anomalous. If the latter is found, i.e. the d-spacing close 
to a transition can be fit to an expression of the form 
[(T — T c )/T c ]~ a , where T c is the critical temperature and a 
the critical exponent, then we intend to quantify the effect 
by determining the critical exponent, and test the theoret- 
ical predictions. Critically behaving bilayer systems are of 
interest because of their low dimensionality and the implicit 
possibility of long range ordering effects in fluid phases. 

The main goal is to study various membrane interactions 
in order to identify the source of the anomalous swelling 
in PC-membranes. We probe the interactions with two 
different sterols, cholesterol and androsten. In terms of 
structure, cholesterol resembles androsten prolonged by an 
extra aliphatic tail. The relative sizes and position of the 
two sterols and the PC bilayer are such that the sterols 
are about half the length of a lipid molecule, and they lie 
towards the water interface 101. 



II. EXPERIMENTAL 
A. Protocol 

Using time-resolved simultaneous small- and wide an- 
gle x-ray diffraction (SAX and WAX), we are able to re- 
trieve structural information in the region about the main 
transition as well as in the neighboring phases. The ex- 
periments were conducted at beam line XI 3 of EMBL 
at DESY, Hamburg. Beam line [|l9| and data acqui- 
sition system [|o| have been been described in detail. 
Here, we report measurements on dimyristoylphosphatidyl- 
choline (DMPC), though we have performed analogous 
experiments on dilauroylphosphatidylcholine (DLPC) and 
dipalmitoylphosphatidylcholine (DPPC). The DMPC and 
cholesterol were from Avanti (Alabaster, AL), and an- 
drosten (androsten-3/3-ol) from Steraloids (Wilton, NH). 
Sample preparation is described elsewhere || . A single ex- 
periment consisted of a first heating from 4 °C to 29 °C , i.e. 
to more than 5 °C above the main transition temperature. 
The sample was then cooled to 10 °C , followed by a sec- 
ond heating up to 29 °C ; i.e. each cycle consisted of three 
ramps. This procedure was chosen to exclude effects from 
improper mixing of lipids and sterols during sample prepa- 
ration, to single out phase effects caused by a slow kinetics 
and/or metastability, and to look at possible hysteretic be- 
havior. The scan rate was 0.5 °C /min within the phases, 
and 0.25 °C /min between 19 °Cand 26 °C covering the 



main transition. This rate was determined in test experi- 
ments with the lipid only. A second sample stayed during 
the experiment in the same sample holder but was not ex- 
posed to the beam. It was always measured as a reference, 
prior to and after the temperature cycle, to check for pos- 
sible radiation damage. 



B. Data 



From our x-ray data we directly determine the positions 
of the first- and second-order lamellar reflections, as well 
as the wide-angle reflections due to the ordering of the 
aliphatic chains. Typical simultaneous small- and wide- 
angle x-ray diffraction patterns of vesicles of DMPC with 
10 mol% androsten are shown in Fig. [p. 
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FIG. 1. Small and wide angle x-ray diffraction patterns 
of dimyristoylphosphatidylcholine vesicles with 10 mol% an- 
drosten. Raw data of the first heating through the main tran- 
sition is shown. The scattering vector s = 1/d — (2sin#)/ A), 
with 2$ the Bragg angle. 

The swelling is seen by the movement of the layer re- 
flections in the SAX-region towards smaller s-values (i.e. 
smaller angles) upon approaching the the main transition. 
The ripple phase is identified by the (lO)-reflection, that 
is seen near the beam stop edge, as well as by the rising 
reflection from the ordering chains in the WAX-region. 

We recorded two orders of reflections in order to 
better determine the position of the reflections since, 
especially in the ripple phase and in the coexistence 
region, there is an overlap due to scattering from 
the in-plane ripple structure. Fig. || shows the 
d-spacing d(T) for the variety of samples measured. 
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FIG. 2. The effects of cholesterol and of androsten on the 
bilayer thickness d of dimyristoylphosphatidylcholine vesicles: 
results of the first heating (see text) are shown. The experi- 
mental points are connected by lines to aid the eye. 

III. RESULTS 
A. Determination of T m from x-ray measurements 

If only SAX data are collected for measuring T m , then for 
PC systems the transition temperature T^ x might be de- 
termined as a major change in slope d d(T)/dT at the low- 
temperature end of the swelling j2| . In particular, when the 
d-spacing of the two phases participating in the transition 
is continuously changing, T m can not be determined with 
precision. In combined time-resolved SAX and WAX stud- 
ies we could define a transition temperature T rn at which 
the intensity of the chain reflection recorded in the ripple 
phase is reduced to half of its maximum. Instead, taking 
into account the coexistence of the ripple phase and the La- 
phase, we here define a transition zone ^.T^ l ax over which 
the intensity drops with increasing temperature. Though 
the WAX measurements describe the chain melting process 
characterizing T m (and the degree of chain ordering could 
serve as order parameter), these measurements depend - 
to a certain degree - on long range periodic order among 
the chains. Moreover, it is by no means obvious that the 
interactions responsible for the chain ordering are directly 
related to the critical bilayer swelling d. One requires, in 
addition, a relation between the cooperativity of chain or- 
dering and critical behavior. 

The d-spacing as well as the intensity of the ripple re- 
flection exhibits a very strong temperature dependence 
near the main transition, as already studied in detail pjl ], 
which is evidence of the non-equilibrium state of the rip- 
ple phase at least near the main transition. The forma- 
tion of the ripple phase itself might thus be a critical 
phenomenon, too. The existence of the ripple phase ap- 
pears to be linked to the processes on the higher tem- 
perature side of the main transition. We also note that 



of all the samples being considered for anomalous behav- 
ior only those displaying a local maximum in d(T) at the 
main transition show a clear (lO)-ripple reflection that 
forms simultaneously during the main transition (Fig. |J). 
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FIG. 3. Contour plot of the intensities recorded in the 
10 mol% androsten (in dimyristoylphosphatidylcholine vesicles) 
experiment as function of the scattering vector s and tempera- 
ture T beginning at 10.4 °C . Small- (first two orders seen) and 
wide-angle (chains) x-ray diffraction patterns are shown. Ti, is 
defined in the text at Table 



B. Is the swelling critical ? 



1. Critical Behavior 



Any singular behavior (power-law- or logarithmic- 
divergent, as well as a finite peak), of a quantity f(T) 
about a critical temperature T c , can be embraced in one 
single functional formula [E2j: 



f(T) ~ - 
a 



(T - T c 



Tr 



- 1 



(1) 



where a is the critical exponent. We preserve the name a 
to point out the similarity of the lattice swelling to the in- 
creasing specific heat (with decreasing temperature) in the 
transition. In a log-log representation of f(T), a fit then 
yields the critical parameters a and T c . If a nonlinear, 
yet well behaved, graph of the log-log data is obtained the 
function f(T) can be described analytically. Therefore, the 
representation of formula (1) serves as a test to distinguish 
between analytical and critical behavior. A critical expo- 
nent of zero needs some further investigation, since it either 
indicates a logarithmic divergence, a cusp-like behavior, or 
a perfectly analytical function with no anomalous behav- 
ior worse than a jump discontinuity p3|| . The log-log-plots 
(Fig. ||) disclose a critical behavior of d(T). 
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FIG. 4. Anomalous swelling: log-log-plot of the bilayer 
thickness d (of vesicles of dimyristoylphosphatidylcholine with 
0, 5 or 10 mol% of androsten) is fitted to Expression (2) . Linear 
fits yield a = 0.01175, 0.02715 and 0.03074 respectively (first 
heating, see text and Table Q). The last point on the line is 
labelled T b . 

Here we address the problem of determining T c for a 
precritical system, i.e. where T c is hidden by the phase 
transition. Formally, this problem is implicit to the defini- 
tion of a critical exponent (Formula 1) which implies the 
knowledge of T c . 



i = {1,6}- The temperature recorded at the data point 
at the low-temperature end (Tb — 23.90 °C) of the thus- 
defined region of anomalous swelling gives the starting tem- 
perature of the T c variation. The system behaves critically 
if all the data points in the anomalous swelling region in- 
deed make a linear function in the log-log-representation; 
if not, the swelling is not critical. In the next step, if the 
data follow a linear function, the slope of this linear func- 
tion (i.e. the critical exponent) is found by minimizing the 
sum of squares of deviations from the line. All of the data 
points enter with the same weight. The first data point 
outside the purely critical regime (at 23.60 °C) is easily 
distinguished, even by visual inspection, by its large devia- 
tion off the linear fit in the log- log representation (Fig. ^). 
The preceding non-deviating point, defining the "break- 
down temperature" Tb, where the critical behavior breaks 
down, is at 23.90 °Cin Fig. ||. The existence of such a 
data point follows from the observation of precritical phe- 
nomena. The procedure finds the global minimum because 
raising or lowering of T c yields a systematic concave or 
convex distortion, respectively, of the string of data points 
from a straight line. 

Table g gives T c , the temperature Tb of the last non- 
deviating data point in the critical regime closest to T c , 
and the critical exponent a, along with the transition tem- 
peratures T™ x obtained from the d(T)-plot, and AT™ ax , 
as well as the amount of anomalous swelling observed 
(d(Tb) — d(28 °C)). These temperatures are illustrated 
in Fig. pi 



2. Determination of T c 

When a precritical system is in the swelling region where 
it displays critical behavior it does not "know" that it will, 
upon approaching T c , undergo a first order transition. So 
we can treat the precritical system over the swelling region 
as a critical system . The fact that the system does not 
"know" a priori the order of the transition it will undergo 
when it swells critically is the base for the theory of critical 
unbinding of membranes Q, which means the d-spacing 
can diverge. Therefore, we can reduce Formula (1) to the 
case a > lEl. 



lnd(T) ~ In 



(T - T c ) 

T r 



d(T) 



(T - T c ) 

T r 



(2) 



In the following, we use Fig. || (DMPC, first heating) as 
an example. To obtain T c we employ a least squares mini- 
mization technique while varying T c . If the system behaves 
critically then all data points (starting from the highest 
temperature data point T) will follow a linear function in 
the log-log plot, upon appropriate choice of T c . The advan- 
tage of this method is that, beginning at the highest tem- 
peratures T (for this example T = 28.05 °C ), it first finds 
the maximum number of data points b that could follow 
a linear function, by comparing the measured d-spacings 
under the constraint d(T) — d(T + i) < d(T + i) — d(T + 2), 
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FIG. 5. Illustration of transition-related temperatures: 
lattice spacing d versus temperature (upper graph, closed 
symbols) measured for dimyristoylphosphatidylcholine vesicles 
with 10 mol% androsten during the first heating. The dashed 
line is a fit to Expression 2. The lower graph shows the maxi- 
mum intensity of the chain reflection (normalized to the inten- 
sity measured at T — 10 °C). The drop in intensity indicates 
the melting of the chain lattice during the main transition. 

The critical behavior is observed (Fig. |i|) in the region 
closest to the critical point over at least a full decade in 
(T — T c )/T c , where data have been collected. For all of 
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the samples listed, the single value of T c itself is to be in- 
terpreted with caution. We remind ourselves that we do 
not directly look at the critical behavior about a second 
order transition point T c . Instead, the first order transi- 
tion intercedes, and thus eliminates from observation the 
region in the immediate vicinity of T c . This region is the 



most important for determining the critical exponent for 
critical phenomena, because there the accompanying long- 
wavelength fluctuations govern the physical properties of a 
critical system and have to be taken into account, e.g. by 
means of RG-theory p3] . 



TABLE I. Critical parameters, T c (critical temperature), Tt (breakdown temperature) and a (exponent), for dimyristoylphos- 
phatidylcholine vesicles with two sterols at various concentrations (cone). Also given are the main transition temperatures from small- 
and wide-angle x-ray diffraction T™ 1 . AT^" 1 , and the anomalous swelling d(Tt) — d(2S °C ) (illustrated in Fig. H). 



Sterol 


Cone. 






T c 


T b 


Tb — T c 


a 


d(T b ) - d(28 °C ) 




mol% 


°c 


°C 


°C 


°C 


°C 




nm 


1 st Heating : 


















None 





23.3 


23.3-23.9 


22.55 


23.9 


1.4 


0.02175 


0.24 


Androsten 


5 


21.2 


21.0-22.3 


21.51 


21.9 


0.4 


0.02715 


0.47 


Androsten 


10 


21.3 


20.5-22.2 


21.47 


21.9 


0.4 


0.03074 


0.54 


Cholesterol 


5 


22.6 




23.03 


23.3 


0.3 


0.01771 


0.39 


Cholesterol 


10 


22.4 


21.9-23.0 


21.24 


22.4 


1.2 


0.03827 


0.42 


1 st Cooling : 


















None 





23.0 


23.6-22.8 


22.45 


23.7 


1.3 


0.02512 


0.27 


Androsten 


5 


21.0 


22.0-19.9 


20.90 


22.1 


1.2 


0.03796 


0.42 


Androsten 


10 


21.0 


21.8-21.0 


21.06 


22.1 


1.0 


0.03409 


0.39 


Cholesterol 


5 


22.2 




21.86 


22.5 


0.6 


0.02714 


0.44 


Cholesterol 


10 


22.1 


22.4-22.0 


19.72 


22.2 


2.5 


0.05164 


0.40 


2 nd Heating : 


















None 





23.2 


23.2-23.7 


22.89 


23.9 


1.0 


0.02264 


0.30 


Androsten 


5 




20.0-22.0 


21.28 


22.1 


0.8 


0.03445 


0.46 


Androsten 


10 


21.3 


21.2-22.3 


21.40 


21.9 


0.5 


0.03329 


0.55 


Cholesterol 


5 


23.0 


22.8-23.4 


22.72 


23.0 


0.3 


0.02453 


0.40 


Cholesterol 


10 


22.3 


21.8-22.9 


20.90 


22.7 


1.8 


0.04173 


0.37 



The advantages of our approach to investigate the crit- 
icality of a system are manifold: no ab initio model needs 
to be assumed, which is a valuable asset especially for crit- 
ical phenomena and phase transitions; no fitting of experi- 
mental data other than the linear fit (to the log-log plots) 
to determine the exponent is needed; and no adjustable 
parameters are used. The breakdown temperature re- 
placing the non-observable critical temperature T c when 
one applies the theory of critical phenomena is not very 
sensitive to the fit. On the other hand, T c and a are more 
sensitive to the fit, also due to the relatively small temper- 
ature range where data have been collected. However, the 
exact value of the critical temperature T c , and hence the 
obtained value of the exponent, is not crucial for the inter- 
pretation of precritical phenomena. We stress that it is the 
breakdown temperature Tb, and not the critical temper- 
ature T c , that is characteristic for precritical phenomena. 
For critical phenomena T b — T c . In our experiments we 
find in all cases Tb > T c . 



3. The precritical parameters 



The exponents derived are all close to zero, and in the 
light of the theory of critical phenomena, a small magni- 
tude of an exponent implies a sharp divergence |^| . The 
exponent a is the quantity that allows us to rank precrit- 
ical phenomena. With one exception a is always smaller 
for the pure lipid than for the sterol samples (Table Q). 
Following the theory for critical phenomena, this implies 
that the d(T) dependence for the pure lipid is the closest 
to criticality, i.e. the pure lipid would display a singular- 
ity of a relatively higher order at T c than would the sterol 
samples. Indeed, looking at the (i(T)-plots obtained from 
the layer reflections (Fig. |^), the pure lipid shows the flat- 
test region between 28 °Cand T™ x (at ca. 23.5 °C). We 
conclude that the higher the concentration of sterol (for 
those samples that display critical behavior) the higher is 
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a, i.e. the less critical is the behavior. This is related to 
the breakdown of criticality further away from T c . 

It is further noticed, for samples showing the strongest 
swelling effect in Fig. | (i.e. d(T b ) - d(28 °C ) for 5 mol% 
cholesterol and 10 mol% androsten), that the two tempera- 
tures obtained in our analysis of the precritical behavior (T c 
and Tb) are relatively close (Table ||). Then, both tempera- 
tures have a qualitatively similar dependence on the sterol 
concentration pointing towards a relationship among each 
other. 



C. Transition-related temperature correlations 

Table | shows that the lowering of T^ x depends on the 
sterol, since androsten lowers T^ x more than cholesterol. 
For both sterols the two concentrations tested, T^ x and 
AT™ ax do not differ significantly However, the number 
of concentrations measured is not enough to conclude that 
these temperatures are only weakly dependent on the sterol 
concentration, as long as the system shows the anomalous 
swelling pH| . There is only a small temperature hysteresis. 
For T c we see a similar behavior with the exception that 
there is an obvious dependence on the cholesterol concen- 
tration. The dependence of T b on the sterol concentration 
is analogous to that of T^ x , but the lowering is not as 
pronounced. AT™ ax features the same tendency as T b . 

The samples displaying the largest swelling effect (i.e. 
d(T b ) - d(28 °C)), which are the 5 mol% cholesterol and 
the 5 and 10 mol% androsten samples (see Table |jj), yield 
a lower T^ x than T c . Of course, this does not mean that 
these samples show a real continuous transition with an 
observable critical point, since we always find T™ x < T bl 
i.e. on cooling, the swelling remains precritical and the sys- 
tem never reaches the maximum swelling before it leaves 
the critical regime. On the other hand the breakdown of 
critical behavior at sub-maximal swelling (e.g. for 5 mol% 
androsten d(T b ) = 6.86 nm< d{T™ x = 6.95 nm) (cf. Fig. 
||) implies that the systems continue to swell below T b but 
not in a critical manner. This non-critical swelling, accom- 
panied by the continued chain ordering (T b is always larger 
than both Tf^ x and the lower end of the transition zone 
AT™ ax (Table is evidence that the chain ordering does 
not obey critical dynamics in the immediate neighborhood 
of the maximum swelling. On heating a decreasing chain 
ordering, combined with a noncritical reduction of the lat- 
tice parameter in the vicinity of its maximum, is seen. T b 
is again always larger than T^ x and the lower end of the 
transition zone AT™ ax (Table ||). Therefore, only at tem- 
peratures T b higher than T^ x does the system move into 
the critical regime (Table ||). For all samples, we generally 
find T^ x on the lower temperature side of the transition 
zone AT™ ax (Table |). Thus, on heating the d-spacing 
does not start to decrease before there is a decrease in 
the chain ordering. On cooling, the swelling does not stop 
before there is already periodic chain ordering established. 
This relationship between the chain ordering and the entire 



swelling - critical and non-critical - seems to be indepen- 
dent of the sterol used and also of its concentration. 



D. Effect of sterols 

1. Membrane swelling 

We see (Table |) that the more pronounced a sample 
shows the swelling effect (i.e. d(T b ) — d(28 °C)) the more 
T c moves towards the high temperature end of the transi- 
tion region AT™ ax . The 10 mol% cholesterol appears to 
not match this rule because the sterol concentration is ob- 
viously high enough to diminish the critical swelling [fyj. 
Only the androsten samples actually have a T c inside the 
transition region AT™ ax . Since the sterols, as basically hy- 
drophobic entities incorporated into the lipid bilayer, inter- 
act in the first place with the aliphatic chains, they decrease 
the tilt angle of the ripple phase ]25| ] upon incorporation. 
As a consequence they also increase the headgroup hydra- 
tion by increasing the headgroup- headgroup spacing |fjT| , 
accompanied by (critical) fluctuations resulting a possible 
softening of the bilayer which in turn could cause an in- 
creased water layer (cf. Model I). One might argue that 
the sterols overpower the origin of the anomalous swelling, 
a PC property, by decreasing the chain and/or headgroup 
tilt. This is not the case since the swelling remains critical 
for smaller amounts of sterols as shown above. A mea- 
sure of the impaired periodic chain ordering (due to the 
presence of the sterols) is the width of the chain reflec- 
tion dependent on the sterol concentration during the main 
transition. The chain reflection widths (measured in detec- 
tor channels) are: 49 (no sterol), 52 (5 mol% androsten), 
58 (10 mol% androsten), 56 (5 mol% cholesterol) and 63 
(10 mol% cholesterol). Because the chains are bound to 
the headgroups, the disorder induced in the chains should 
be transferred to the headgroups as well. Apparently, the 
impaired chain ordering, due to the sterols, yields a larger 
exponent a indicating a weakened criticality (Table |) but 
is not a monotonous function of the sterol concentration. In 
other words, the impaired periodic chain ordering (width) 
mirrors the degree of criticality a, but is not the origin of 
the anomalous swelling. 



2. Cholesterol versus androsten 

For the cholesterol samples our data do not confirm the 
mutual exclusion of critical behavior and chain ordering 
because T b falls into the transition range AT™ ax . How- 
ever, the tentative statement made above - that the chain 
ordering is not the origin of the critical behavior - remains 
valid, last but not least because increased chain ordering 
does break the criticality in the cholesterol samples, too, as 
evidenced by a T b within the transition range AT™ ax (cf. 
Table ||). Thus, swelling in the presence of cholesterol can 
exist with limited periodicity in chain ordering. Possible 
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explanations for the differing amount of chain ordering, 
required to break the critical behavior observed between 
the two sterols are (a) a sterol-specific interaction with the 
lipid molecules, and (b) an active domain-boundary seek- 
ing of sterols in the coexistence region. Sterols, as sensors, 
have revealed an imperfect coupling of the chain order- 
ing processes with the layer dynamics, i.e. the swelling, 
pointing towards a phenomenon affecting the hydrophilic- 
hydrophobic interface. Since none of the four models ex- 
cludes changes in the headgroup-hcadgroup interaction, we 
now discuss the effects of critical fluctuations and hydra- 
tion on these interactions. 



IV. DISCUSSION 



A. Testing the four Models 



In the remainder of this paper, we associate the experi- 
ments with the four models discussed earlier 0. Because 
undulations and fluctuations are basic to Model I, we first 
discuss critical fluctuations and the correlation of fluctu- 
ations with swelling. The effect of chain ordering on the 
breakdown of criticality is connected with Model II, and hy- 
dration is connected with Models III and IV. The advent 
of a ripple phase appears to be connected with the main 
transition and pre-criticality. Finally, thermodynamic and 
nonequilibrium effects are discussed, and the importance 
of low scan rates is stressed. 



B. Critical fluctuations 



1. Fluctuations in membranes 



Lipowsky and Leibler |1J] predict, on the basis of Hel- 
frich undulations, that the bilayer thickness d should be lin- 
ear in 1/(T-T C ) (Model I), i.e. a = 1. Our results (Table 
|lj) give exponents a that differ significantly from 1 , whereas 
earlier experiments on lipid membranes |6|] and other soft 
condensed matter systems Q supported the predictions 
|T^| . The magnitude of our precritical exponent a is rea- 
sonable, since an exponent of the order of 1 would yield 
a much larger swelling effect over the entire experimental 
critical regime of 28 °C down to T^. Experimentally, even 
the largest amounts of critical swelling observed (Table |, 
d(Tb) — d(28 °C )) do not reach a tenth of the unit cell size d 
(of around 7 nm). When plotting our data (Fig. ||) against 
1/(T — T c ), we observe clear deviations from the predicted 
linear behavior a single linear fit is not reasonable. 




FIG. 6. Critical bilayer thickness (of vesicles of 

dimyristoylphosphatidylcholine, first cooling, scan rate 
0.1 °C/min, T c = 22.55 ± 0.05 °C); a prediction [TJ gives 
a straight line. 

Thermal fluctuations trigger the main transition. To in- 
vestigate the role of enhanced fluctuations (all Models) on 
the swelling we consider the transition-related temperature 
derived from the layer reflections T™ x , since it marks the 
temperature at the end of the swelling, critical plus non- 
critical. Our data (Table ||) show that the addition of sterol 
shifts this transition temperature T^ x . This means, only 
a partial periodic chain ordering is required to obtain the 
maximum d-spacing in the transition region. Even more, 
for the samples showing the strongest swelling effect yield- 
ing a local maximum in d(T) (Fig. ^, 5 and 10 mol% 
androsten), the continued ordering upon lowering temper- 
ature reduces d{T) slightly. A more pronounced local max- 
imum in the transition range is seen in the pure DMPC as 
well, but here it covers a broader temperature range com- 
pared to the two sterol samples. This peculiar relationship 
between chain ordering and d-spacing leading to the local 
maximum might be explained by the changing contribu- 
tions of in-plane and out-of-plane fluctuations. Whereas 
the in-plane fluctuations affect the periodic ordering of the 
chains, the out-of-planc-fluctuations in first place deterio- 
rate the periodicity of the layer order. The lower energetic 
in-plane fluctuations are expected to contribute more at 
the lower temperature end of the swelling when the inter- 
layer spacing becomes larger and larger and neighboring 
membranes are basically decoupled. However, when the 
chains start to order due to increasing chain-chain-coupling 
these fluctuations will then be suppressed. In the theory of 
liquid crystals -lipid model membranes can be considered 
as such, in-plane fluctuations affect the bending modulus 
K c and the out-of-plane fluctuations affect the compression 
modulus B [[La]. Here, fluctuations broaden the width of 
the x-ray reflections as a function of the order caused by 
the decreasing correlation length of the scattering units. 
Fluctuations, or defects of the first kind, also reduce the 
maximum intensity with increasing diffraction order pTj . 
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2. Fluctuations: Correlation with swelling 



C. Effect of chain ordering on the breakdown 



The intensities of the x-ray reflections give information 
about membrane fluctuations. Our data show (Fig. 0(a)) 
a decrease in the ratio of the maximum intensity of the 
first order to the second order reflections, with increasing 
sterol concentration and with decreasing temperature, in 
the swelling region (of 28 °Cdown to TJ,). 

Similarly, Fig. ^(b) shows that the ratios of the respec- 
tive integral widths (integrated intensity/maximum inten- 
sity) has approximately the same variation in T for all 
sterol concentrations for T > 
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FIG. 7. Ratio of maximum intensities (a) and integral widths 
(b) of first and second order reflections as function of tempera- 
ture T (Tb: breakdown temperature) for vesicles containing up 
to 10 mol% sterol; symbols as in Fig. ^. 



When comparing the breakdown temperature TJ, to the 
main transition zone AT^ ax (Table |) we observe that 
virtually coincides with the upper bound of AT™ ax for the 
pure lipid and the androsten samples on heating and cool- 
ing. For the cholesterol samples this is clearly not so (Table 
|): A certain degree of chain ordering appears to be per- 
mitted as Tb moves into the transition zone AT™ ax . This 
ordering effect must be related to the difference between 
the two sterols and their interaction with DMPC. The fact 
that for the pure lipid Tb coincides with the upper bond 
of periodic chain ordering (given by AT™ ax ) is persuasive 
evidence for the following scenario. On the one hand, Tb 
marks (by definition) the observed lower temperature limit 
of the critical behavior (as defined by Formula 2) in the 
L Q -phase. On the other hand, the upper bound of AT™ ax 
denotes the higher temperature limit of the periodic chain 
ordering of the ripple-P^'-phase. Hence, on cooling, the 
onset of chain ordering breaks criticality (as is illustrated 
in Fig. H where the Tb line intersects with the edge of a 
detectable WAX reflection.) Therefore, changes in inter- 
actions within the hydrophobic membrane core, i.e. the 
steric repulsion and the van der Waals interaction between 
them, coupled to the chain ordering can not directly control 
the critical swelling, so Model II is not supported. (Mod- 
els I, III and IV are indifferent towards chain ordering ef- 
fects.) Neither the steric repulsion nor the van der Waals 
interaction among the chains are expected to change dras- 
tically when entering the coexistence region, because both 
do not depend on the degree of periodic order. This is sup- 
ported by the fact that one only observes critical swelling in 
PC but not e.g. in phosphatidylethanolamine (PE) mem- 
branes. (Additional data taken on mono- and dimethylated 
PE do show the swelling.) 



In light of the theory of x-ray diffraction, the similarity 
of the dependence of both integral width and maximum 
intensity ratios on sterol concentration is to be interpreted 
as an increase of static defects p7|] , whereas the decreasing 
ratio of the maximum intensities with decreasing temper- 
ature when approaching Tb (i.e. with increasing swelling) 
indicates a reduction of fluctuations (akin to a reduction in 
the Debye- Waller factor). Since the reduction is relatively 
stronger in the sterol samples than in the pure DMPC, 
it follows that sterols reduce fluctuations. In consequence 
we have to exclude a significant fluctuation-driven contri- 
bution to the critical swelling, i.e. Model I is not sup- 
ported. (However, they increase the criticality of the sys- 
tem as expressed by a (Table |)). The observed fluctuations 
are therefore not critical fluctuations and can be described 
in terms of the classical Ornstein-Zernicke theory ^8|; in 
other words, the system is sufficiently far away from a crit- 
ical point that a Landau-type theory is appropriate for the 
description of the fluctuations. 



D. Hydration 

The magnitude of the critical swelling d(Tb)~ d(28 °C ) is 
increased in the samples containing a small amount of sterol 
(< 10 mol%) (Fig. I and Table |), compared to that of the 
pure lipid. The adding of small amounts of sterol does 
not significantly decrease the area per headgroup (^9|, so 
that the difference in area per molecule in the neighboring 
phases remains large (about 15% -20% §). In tur 11 we ex- 
pect an increased amount of water in the headgroup region 
possibly leading to a change in headgroup conformation, as 
is proposed by Model IV. NMR measurements indicate that 
water molecules might penetrate past the hydrophilic head- 
group towards the hydrophobic membrane core Jl6| . The 
changing hydration affects the hydrophilic/hydrophobic in- 
terface. The dynamics of an interface is known to depend 
on the local curvature of the interface |30| and can be- 
have critically under certain conditions [pl|. The addition 
of 15% sterol does not provide appropriate conditions to 
observe critical swelling: Due to phase separation effects 
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expected (a miscibility gap Q] and the relatively low tem- 
perature in these experiments, no critical behavior can de- 
velop (cf. Fig. ||). Also, the diffraction technique is insen- 
sitive to short-range order. We have studies in progress on 
the miscibility gap p2fl . A change of the hydration shell 
about the headgroup is thus a candidate for the critical 
swelling effect. We further remark that a decreased head- 
group area upon adding sterol does increase the amount of 
swelling but reduces the criticality of the swelling (Table 
If). One is tempted to identify the kink in our results (Fig. 
|), at around 24.55 °C, as the boundary of a critical hy- 
dration dominance. (The changes in slope in Fig. ^| may 
be interpreted as a critical slowing down ]3q] .) Lipowsky 
and Leibler imposed a hard confinement scenario in their 
derivation jl4| , which is perhaps not valid here |34| , where 
a change in environment can be attributed to a change in 
hydration. Also, NMR data jl6| indicate a change in hy- 
dration at about this temperature. Moreover, dynamical 
scaling can break down despite the persistence of interfa- 
cial undulations j35| . To clarify the difference between the 
critical behavior of a changing headgroup hydration (Model 
IV) and bilayer thickening (Model II), we point out that the 
hydration depends on the headgroup environment whereas 
the bilayer thickening in Model II is only a property of the 
chain conformation. The headgroup environment and the 
hydrophobic chains are coupled, and thus headgroup and 
hydrocarbon chains could be considered as an integral en- 
tity (Model II) . The differences between the two models are 
(a) the coupling is not stiff, so the behavior of chains and 
heads can be different and might not be mapped directly 
onto each other, and (b) chain-chain interactions involve 
lipids only, whereas hydration dynamics requires solvent 
(here water) and lipid molecules. These differences could 
give a clue why a lipid system can possess only precritical 
behavior, but not critical behavior. 

When entering the coexistence region AT™ ax the head- 
groups feel the ordering of the coupled chains. As shown 
above the chain ordering correlates with the breakdown of 
the critical swelling and hence breaks a critical interface dy- 
namics. In response the system finds a new (metastable) 
equilibrium position forming the ripple structure. It has 
been suggested that the swelling of the hydration shell 
at a hydrophobic/hydrophilic interface is a quite general 
phenomenon f36| . Also, osmotic pressure experiments on 
various PCs seem to show a critical dependence of the d- 
spacing on the osmotic pressure |34| . 



V. SUMMARY AND CONCLUSION 

We demonstrated that the anomalous swelling in phos- 
phatidylcholine membranes, of 14 carbons chain length, is 
precritical. By analyzing the nonlinear d(T) behavior - 
in terms of the theory of critical phenomena - observed 
for the pure lipid and various sterol concentrations in the 
neighborhood of T c , we were able to link the anomalous 
swelling to a changing headgroup hydration. 



We propose that the formation of the ripple phase is 
the counterpart on the low-temperature side of the tran- 
sition to the anomalous swelling on the high-temperature 
side. It is suggested that the ripple phase and the anoma- 
lous swelling are caused by the hydration dynamics of the 
headgroups. 

We have found no evidence of critical entropic bilayer 
undulations quoted in Model I. Model II suggests a critical 
thickening of the bilayer. The increase in the conforma- 
tional ordering of the chains in the fluid phase, as they 
approach the coexistence region, is not known to obey crit- 
ical dynamics. 

Our sterol experiments show that adding small amounts 
of sterol does not break - but reduces - the criticality of 
the system, which seems to be linked to the larger head- 
group area at the interface. This seems to be in agreement 
to Model III, suggesting that a decreased interfacial surface 
increases the hydration force j37| , however a decrease in in- 
terfacial area over the region of anomalous swelling can not 
be supported by our data. This suggests that the hydra- 
tion force does not change in a critical manner. The mech- 
anism proposed by Model IV, a conformational change of 
the headgroups due to changing hydration, is supported 
by our experiments and can explain the critical behavior. 
However, the model does not explain the amount of swelling 

From our data the following picture emerges: the chang- 
ing hydration of the headgroups is responsible for the criti- 
cal behavior and indirectly coupled chain ordering dynam- 
ics provides the full magnitude of the swelling. Further, 
the similarity of the transition-related temperatures deter- 
mined from the SAX- and WAX-reflections (Table |) are 
evidence for a relationship of the conformational chain dy- 
namics to the systems' transition from the fluid to the rip- 
ple phase, but not for the transition from the critical to the 
coexistence region, which is found to be related to the long 
range periodic chain ordering dynamics. 

Whereas DMPCholine exhibits critical swelling, DM- 
PEthanolamine does not |58|. The much larger polar 
choline headgroup, the only structural difference between 
both lipids, is easy to hydrate. The existence of lattice 
critical swelling may well be connected with this differ- 
ence in changing the hydration, and with the differences 
observed, between DMPC and DMPE, in differential scan- 
ning calorimetry [ |39"| of the main transition. 

We emphasize that our precritical exponents describe the 
swelling behavior over the entire temperature range (where 
swelling is observed, in all of our samples) supporting the 
subtle interface dynamics picture. However, a necessary 
requirement for an ab initio theoretical description is the 
inclusion of the formation of the ripple phase. This work 
suggests that the critical behavior is caused by the hydra- 
tion dynamics. 
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